o 

(N 



Mon. Not. R. Astron. Soc. 000, [THS] (2010) Printed 7 July 2010 (MN M^X style file v2.2) 

The influence of modified gravitational fields on motions of 
Keplerian objects at the far-edge of the Solar System 

Michael Sokaliwska 1 *, Hans-Jorg Fahr 1 ! and Pavel Kroupa 1 ! 

1 Argelander Institut fur Astronomie, Auf dem Hiigel 71, Bonn, Germany 
MNRAS, accepted 



in 

w 
6 



> 

in 

q 

o 
o 



X 



ABSTRACT 

We investigated the impact of three different modifications of Newtonian gravity on 
motions of Keplerian objects within the Solar System. These objects are located at 
distances of the order of the distance to the Oort cloud. With these three modifications 
we took into account a heliocentric Dark-Matter halo as was indicated by Diemand 
et al., Modified Newtonian Dynamics (MOND) and a vacuum- induced force due to 
a locally negative cosmological constant A_ derived by Fahr and Siewert. In gravi- 
tationally bound systems it turns out that all three modifications deliver the same 
qualitative results: Initially circular orbits for the pure Newtonian case are forced to 
convert into ellipses with perihelion migrations. The quantitative consideration, how- 
ever, of the orbital parameters showed strong differences between MOND on the one 
side, and Dark-Matter and A_ effects on the other side. 

Key words: Physical Data and Processes: Astrometry and celestial mechanics: grav- 
itation, celestial mechanics, Solar System: Oort cloud, methods: numerical 



1 INTRODUCTION 

i In modern astronomy there are some enigmas that physi- 
cists have not been able to solve yet. These enigmas 
seem to indicate the necessisty to depart from Newtonian 
gravity. One of th ese enigmatic phenom ena is the so called 
Pioneer- Anomaly (jAnderson et al.l 12002 ) that seems to be 

i an anomalous acceleration of the two space-craft, Pioneer 
10 and 11, that were launched in the 1970's and move 
into two opposite directions of the Solar System. These 
spaceprobes show an acceleration towards the Sun of the 
order of 10~ 10 msec" 1 that seems to be independent of 
Solar distance, at least over Heliocentric distances between 
20 and 70 AU, and cannot be explained by usual Newtonion 
dynamics. It seems that there exists an additional force 
responsible for this phenomenon, like modified gravity, 
but the constancy of the acceleration poses an important 
challenge. 

The second indication for modifed gravity manifests 
in the galactic rotation curve problem. The observed flat 
rotation curves at large galactic radii are not explainable 
by Newtonian dynamics alone. 

This paper is inspired by these phenomena. We present 
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in Section 2 of this paper three different modifications of 
gravity. Section 4 deals with the results of numerical calcu- 
lations of the impact of these modifications and illustrates 
how these modifications act. The last section is a short out- 
look and suggests further possible investigations. 



2 MODIFICATIONS OF GRAVITY 

Because Dark-Matter (ACDM) has become an essential in- 
gredient of the standard scenario in modern cosmology the 
first modification that is investigated in the context of this 
paper is a helio c entric Dark-Matter halo as was indicated by 
iDiemand et al.l 1J2005T ) from Neutralino Dark-Matter simula- 
tions. 



2.1 Heliocentric Dark-Matter minihalo 

In order to investigate the influence of a halo component of 
Dark-Matter it was necessary to choose a density-profile for 
the halo that was assumed to be bound to the Solar System 
and centred at the position of the Sun. 



According to lAnderson et al.l (120021 ) who investigated 
whether a DM halo could be responsible for the Pioneer- 
anomaly we decided to use a halo that causes an additional 
acceleration towards the centre that is constant and indepen- 
dent on heliocentric distances, in view to indications given 
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Figure 1. Deviation in radial distance r for different Heliocentric 
DM halo-masses: One orbit of Uranus (« 84 yr) has been inte- 
grated for both cases, for the usual Keplerian case and the case 
taking a --Dark-Matter density profile into account. The differ- 
ences in the absolute distance between the two cases has been 
calculated for every time-step. The enclosed mass is given for a 
radius of 50 AU 



by the Pioneer-anomaly. We started with a halo-density- 
profile pDM oc r" where n is a real number, 



Pdm{t) = par"" 



(1) 



Substituting this expression into the spherical Poisson- 
equation 



L d ( 2 d 
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dr 
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(2) 



leads to the following expression for the gravitational force 
connected with the DM-halo 
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n + 3 

If this force is constant over heliocentric distances, n has to 
be chosen equal to — 1 and the integration constant C\ — 0. 
This then leads to the following density-profile 

rppo 
r 



p(r) 



(4) 



where ro — 1 AU and po is the density at ro, that causes a 
constant additional acceleration towards the centre. 



2.1.1 Halo-mass 

The question is how the halo-mass can be selected such 
that the orbits of the outer planets, l ike Uranus, persist? 
For that purpose I Anderson et alj 112002 ) used data from the 
ephemeries of DE200 ( Anderson et al.l ll995T ) that contain 
Uranus' ephemeries down to an accuracy of 10 -5 AU in ra- 
dial distance. The halo-mass enclosed within a given dis- 
tance can be calculated via 



Mdm —4n r p DM {r)dr. 



(•») 



With that the halo-mass is fixed by the value of po- Cal- 
culating the deviation of the orbit of Uranus from the pure 
Keplerian case for different po leads to the results shown in 
Fig. 1. 



It is obvious that a radial distance deviation of the order 
of 10 AU is given for an enclosed halo-mass, within 50 AU, 
< 10 -6 M Q that corresponds to p « 4.4 x 10 5 M pc -3 . 
Such a halo-mass is in go od agreement wi t h the result s from 
an investigation done by I Anderson et alj ([1998, 2002). 

Using this enclosed mass delivers a radial acceleration- 
component of the order of 10 -12 msec -2 which is two 
orders of magnitude smaller than the measured anomalous 
Pioneer-acceleration. 

Parenthesis: Comparison with a standard-model 
Is a Dark-Matter halo with a density-profile as given by eqn. 
(4) in good agr eement with a genera l Navarro- Frenk- White 
(NFW ) profile l|Diemand et alj|2005l ; lHernquist|[l990l ; IZhaol 
Il996h ? The NFW-profile is given by 



p(r) = 



pO,NFW 



i+ter 



(?—/)<* ■ 



(6) 



where r s is the scale radius of the halo and po,nfw is a 
constant density. For a = 7 = 1 and /3 = 3 this profile 
reduces to a standard NFW-profile 



p(r) = 



P0,NFW 



(7) 



In the limit of small distances compared to the scale 
radius (r << r s ) this profile leads to the desired — 
profile. Compared to Pioneer-like distances (20 — 70 AU) 
r s has to be choosen very large, e.g. for r — 50 AU, 
r s > 5000 AU. With this the constant accelera- 
tion a — 2irGr s po,NFW ~ 10 -12 msec -2 leads to 
Po,nfw < 47 M pc -3 . 

Using a correlation between realistic Halo-parameters 
(Darbringhausen et. al., in prepration), a halo with a scale 
radius r s > 5000 AU would have po,nfw < 30 Mqpc" 3 . 
This value and the one stimated for a scale radius larger than 
5000 AU are of the same order of magnitude and it seems 
that the used density profile (4) is a realistic approximation. 



2.1.2 DM Critical Radius 

The distance from the centre where the Newtonian gravi- 
tional force due to the central mass (the Sun) equals the 
radial acceleration due to the halo with a density-profile 
according to (4) is called the critical radius. This radius 
turns out to be of the order of r cr tt ~ 50000 AU, much 
larger than distances of about 70 AU, so the gravitational 
contribution of such a halo to the central gravity is negligi- 
ble at such distances and would not lead to large daviations 
from pure Keplerian motion. 



Objects expected to move at distances of the order of 
this critical radius are located in the Oort cloud. This cloud 
is assumed to be a reservoir of comets that represents the 
outer most material border of our Solar System and lies at 
a distance of about 50000 AU. 

So it is convenient to consider the influence of the 
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DM-modification on Oort-cloud objects. The calculated ad- 
ditional radial acceleration was included into the used orbit 
integrating program that will be presented in Section 3. 



2.2 Modified Newtonian Dynamics (MOND) 

2.2.1 Newtons second law 

iMilgroml |l983j) introduced a function fj,(z) in Newtons 
second law where z — — , with the 'actual' acceleration a 
devided by some critical acceleration ao which is a universal 
constant fa 2 x 10~ msec - , derived from galactic rotation 
curves and turns out to be of the same order of magnitude as 
the Pioneer-acceleration. With this function Newton's sec- 
ond law is written 



/ a\ 
F = man, I — I = mojv 
Veto/ 



(8) 



where a is the actual acceleration of the object while 

GMp, 
9N = 



is the strict Newtonian acceleration. 



The function n(z) is not further specified but has the 
limiting-values 



fi(z) = 



if z > 1, 
if z < 1, 



(9) 



The second case corresponds to an actual acceleration a that 
is much smaller than the critical acceleration ao. This leads 
to an expression for a 



mafj, 



( d \ a<^a 

\a J 



such that 
a = y/gNdo. 



I'd 



— mgN, 



(10) 



(11) 



This expression does not describe a continuous transi- 
tion from the pure Newtonian regime (a >> ao) to the pure 
MONDian regime (a << ao). So we chose a function that 
fullfills the required limit, 



fi(z) = 



1+2 



With this function equation (8) leads to 



9n , t g% , 

— + I/ — +9NCL0, 



(12) 



(13) 



which is an expression for the actual acceleration that de- 
scribes a continuous transition between the two regimes. 
This was the second type of acceleration considered for com- 
parison studies in the orbit-integration routine. 



2.2.2 MOND-Critical radius 

In anology to the DM-case we calculated the radius were 
the gravitational acceleration due to the central mass equals 
the criticial (characteristic MONDian) acceleration ao to be 
Vera ~ 5000 AU. So also in the MONDian case there is no 
observable influence on Keplerian motion expected on scales 
of the planetary Solar System (sa 50 AU). 



2.3 Vacuum-type force: Similar to a negative 
cosmological constant 

The third modification of gravity considered here is a 
vacuum-induced force connected with a 'negative local 
cosmological constant' A_. 

Our Universe is assumed to be describable by the 
Robertson- Walker- metric 



a (t)(dr +r (dr + sin 6dcj> )), 



(14) 



where a is the scale parameter of the universe, c the velocity 
of light, r the distance and dt an infinitisimal time-interval. 
6 and <j> denote the polar and the azimuth-angles. This 
metric leads to an exact solution of the Einstein Field 
equations and describes a homogeneously matter-filled and 
isotropically curved universe expansion or contraction. 

However, gravitationally bound systems that have 
been decoupled from the homologous cosmic expansion due 
to gravitational collapse should show a different metrical 
dynamic behaviour than the residual universe. 



iFahr fc Siewerj 1)2008! ) picked up this idea in order to 
describe the evolution of gravitionally bound systems like 
our Solar System, and started with an equation of motion 
given by 



d 2 x 
dt 3 " 



GM a(t) 
' X 2 + a{tf' 



(15) 



where x describes the distance from the centre of the Solar 
System, G is Newtons gravitational constant and M the 
effective central gravitating mass. The first term on the 
RHS results from the attractive gravitational field whereas 
the second term corresponds to a repulsive force connected 
with the global expansion of the surrounding u n iverse . 
This equation was derived by ICooperstock et al.l (| 19981 ) 
and describes the impact of the cosmic expansion on the 
dynamics within gravitationally bound systems (the 'local 
systems'). 

Describing the extent of a local system by a "local scale 
factor" l(t), that corresponds to the scale factor of a local 
RW-metric and considering the above equation at the border 
of this system leads to 



GM ,,'a(t) 



(16) 



This second order differential equation can be trans- 
formed into a first-order onqj by multiplying with 21 which 
leads to 



d ;2 , , „ d GM a d , 2 

— fit) = 2——— h -r. 

dt v ' dt I adt 



(17) 



It ca n formally be in t egrate d using a relation also de- 
rived by IFahr fe Siewertl (|2008l . see Appendix there) and 
leads to 



P = 2GM (r 1 



rM + 

b rec I i 



Wa + a_\ (/2 
2 V a a rec J y 



I 2 )+i 2 

''Tec I ' "re 



.(18) 



1 The co mplete derivation of th e following formalism can be re- 
viewed in lFahr fe Siewed ||2008|1 . 
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where l rec denotes the local scale factor at 
recombination-time t rec and a rec the global scale fac- 
tor at this time. 

Solving this equation in its general form needs numerical 
methods and is not trivial. So it is convenient to analyze the 
asymptotic behaviour of this equation. At first consider the 
case that the scale factor I at time t is much larger than the 
associated scale factor at recombination time, l(t) S> l,-ec- 
Then the above equation reduces to 

P = _^ + 1(« + M P + L . (19) 

The RHS of this equation must always be positive and leads 
to the condition 



u, (lrec 



2GM 1 fd + Orec^ ,2 ,2 

Lrec « 



I 4- / 



(20) 



This condition then suggests that, if fullfilled, a region of 
'"gravitationally"' bound spacetime may under such condi- 
tions only grow large, if the central mass is sufficiently small. 

In the opposite case, l(t) <C lrec, the equation of motion 
of the local scale factor yields 



I 1 



2GM l/o oV 



I 



1 /S drec \ ,2 p 

\a a rec / 



which leads to the condition 



2GM 



+ 1 > 



1 /a a 



/ (2 ftrec \ 

Va a rec / 



(21) 



(22) 



and can be interpreted as meaning that a sufficiently large 
central mass causes a n asymptotically small v acuole of the 
Einstein-Straus type IjEinstein fc Strauslll945r ). 

Using solutions of the well-known Friedmann-equations 

t \ 2 / 3 

(23) 



a(t) — a re c ( - — J 

equation (18) leads to an equation 






with the parameter functions 
a = 2GM, 

7 (t) = 



(24) 



<&lrec \ ZLtJVl Ipec ( -i ^rt 



! + ■ 



f 2 



(25) 



1 / 1 + « 



f 2 



that formally is very similar to the conventional Friedmann- 
equation that is given by 



CD 



a\- H$n K A£_ 



(26) 



The comparison with the conventional Friedmann- 
equation shows a pronounced relationship between equation 
(24) and (26), despite some important differences between 
these two equations. 

The first term on the RHS of Eqn. (24) corresponds to 
the cosmological mass density term; the second term with 



P(t) seems to correspond to the curvature term with the dif- 
ference that this curvature is time-dependent here. The last 
term corresponds to a vacuum energy term, but is negative 
and time-dependent. This special result will be investigated. 



Assuming the relation 7(f) 



leads to an equiv- 



alent of the "cosmological constant" that is today (t 



13, 7 x 10 9 ys and t re 
44 1 



3 x 10 5 yr) 



A 



-1 x 10" 



(27) 



As showed in I Westermannl (|2004r ) the force that is act- 
ing on a particle, taking the influence of the vacuum into 
account, is given by 



GM c 2 A 
H 3 

or in the case here, 

GM c 2 A_ 



+ 



(28) 



(29) 



r- 3 " 

where A_ is the above derived negative constant which leads 
to an additional attractive force towards the centre. From 

this one can derive the corresponding potential via J^r- = 

F(r), 

,, , GM c 2 A„ 2 ,o„n 

$(r) = 1 r + const. (30) 

r 6 

The attraction derived fro m this potent i al inc reases 
with increasing distance oc r 2 . IWestermannl (J2004) con- 
cluded that there is no observable influence of A expected 
on scales of the inner Solar System. However, the above de- 
rived NEGATIVE cosmological constant A_ is about three 
orders of magnitude larger than the conventional A and thus 
changes earlier conclusions. 



2.3.1 A_ Critical radius 

The force due to the negative vacuum energy density is 

c 2 A_ 

— j= — r. Since A_ is negative the induced force is directed 
towards the centre, just as the gravitational force. 

The above expression leads to the condition for the crit- 
ical radius, 



_ / 3GA/ V 

" \c 2 A- ) 



(31) 



It turns out that this radius is about r cr u ~ 4.5 x 10 AU. 
At this distance from the Sun the influence of the vacuum 
cannot be neglected anymore, or, in other words, it is not 
possible to use perturbation theory because the "perturba- 
tion" is neither small nor time limited on a small interval. 
By comparison, the general cosmological constant with an 
upper limit of 10~ 47 m -2 would lead to a critical radius of 
about 8 light year. 

Obviously a distance of the order of the critical radius 
is too large to expect a direct influence on planetary orbits. 
But if one takes the Oort cloud into account with an extent 
of about 1 light year, it should be possible to see deviations 
from pure Keplerian expectations. 

Since the critical radii of the three gravitational modifica- 
tions that are investigated in this paper, are very large, they 
are not able to explain phenomena like the Pioneer-anomaly. 
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Figure 2. Resulting orbits for a distance of 50.000 AU under the 
influence of a Dark-Matter minihalo. The integration time was 
14.02 MYr. 
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Figure 3. Resulting orbits for a distance of 28.000 AU under the 
influence of a Dark-Matter minihalo. The integration time was 
12.65 MYr. 



3 THE ORBIT-INTEGRATING PROGRAM 

In order to investigate the impact of the three modifications 
of gravity derived above we use an orbit-integrator based 
on an embedded Runge-Kutta method of 4th order. The 
energy and angular momentum turns out to be a conserved 
quantity for all modifications. 

Since the distances that were considered here are of the 
order the distance to the Oort cloud (« 10 4 AU), the whole 
problem reduced to a two-body problem: The central mass 
and the orbiting test-object. 

We considered eleven cases with initial conditions that 
would cause clear circular orbits, -with initial distances in 
a range between 5000 — 50000 AU- in a pure Keplerian 
potential and two cases that correspond to ellipse^]. 

The integration times were chosen differently (see Fig. 
2 - 7) in a range between 4.9 and 17.5 MYr. 

It should be noted here that during the investigation of 
these three modifications every influence of external fields, 
like the Galactic tidal field, has been neglected. So only the 
gravitational force and the additional forces have been con- 
sidered. 



4 RESULTS OF THE ORBIT-INTEGRATION 

Fig. 2-7 show the results of the orbit-integration for all three 
modifications. The solid line marks the circular orbit in a 
pure central-mass Keplerian case and the dashed line de- 
notes the orbit under the respective gravitational modifica- 
tion. 

As one can see all three modifications cause strong devi- 
ations from pure Keplerian orbits around the central mass. 
Circular orbits become elliptical ones that in addition expe- 
rience an aphelion-migration, a so called rosette-orbit. This 
qualitative behaviour occurs in all cases. The differences 
lie in parameters like the value of the aphelion-migration 
and the eccentricity for one turn. The results of the investi- 



2 These cases showed the same qualitative behaviour as the cir- 
cular cases. 
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Figure 4. Resulting orbits for a distance of 50.000 AU under the 
influence of MOND. The integration time was 12.3 MYr. 



gation of these parameters are presented in the next section. 



4.1 Orbital parameters 

In order to compare the impact of Dark-Matter, MOND 
and A_ quantitatively we introduced two kinds of orbital 
parameters, the 'passive' orbital parameters and the 'active' 
ones. 




" 3 %0 -20-10 10 20 30 
x[kAU] 

Figure 5. Resulting orbits for a distance of 28.000 AU under the 
influence of MOND. The integration time was 4.9 MYr. 
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Figure 6. Resulting orbits for a distance of 50.000 AU under the 
influence of A_. The integration time was 17.5 MYr. 
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Figure 7. Resulting orbits for a distance of 28.000 AU under the 
influence of A_. The integration time was 13.6 MYr. 



The passive orbital parameters describe the orbits the 
test-object would be thought to move on under the assump- 
tion that no additional force or modification of gravity is 
valid. An observer who assumes a pure Keplerian potential 
would identify these parameters with the observed spatial 
and velocity components of the object. These parameters 
were calculated at every time step giving the true position 
and velocity values of the object. The passive eccentricity 
and the semi-major distance were the two passive orbital 
parameters considered here. 

The second type of orbital parameters, the active ones, 
describe the actual shape of the orbit caused due to the 
gravitational modification. So the active eccentricity is the 
one of a fitted reference ellipse approximating one turn of 
the rosette-orbit, whereas the aphelion-migration describes 
how much the orbit turns per day. 



4-1.1 Passive orbital parameters 

The used orbit-integrating program was able to calculate in- 
stantaneously the passive orbital parameters for every time- 
step. The passive semi- major distance is given by 

(32) 



2 — rv 2 /v' 

where a is then the passive semi-major-distance, r the dis- 
tance from the centre, v the velocity and v = GMq Newtons 
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Figure 8. The behaviour of the instantaneous best-fitting (pas- 
sive) semi-major-axis (SMA) and the aphelion distance with time 
under the influence of a Heliocentric Dark-Matter mini-halo. For 
the initial conditions a Newtonian circular orbit with 50000 AU 
radius was used. 
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Figure 9. The behaviour of the instantaneous best-fitting (pas- 
sive) semi-major-axis (SMA) and the aphelion distance with time 
under the influence of MOND. For the initial conditions a New- 
tonian circular orbit with 50000 AU radius was used. 



gravitational constant times the central mass. 

The passive eccentricity was then calculated via 



e=\/l-— , 
av 



(33) 



where p^ denotes the constant angular momentum. 

The calculation for the passive orbital parameters was 
done for one orbit with the initial conditions in velocity- 
and space-components corresponding to a circular Newto- 
nian orbit with 50000 AU radius. Fig. 8-11 show the time- 
dependent behaviour of the semi-major-distance for all three 
modifications. 

As one can see in all three modifications the semi-major- 
distance oscillates. In the case of Dark-Matter and A_ the 
values of a stay positive for all times. Not so in the MON- 
Dian case. 

In MOND the semi-major-axis 'jumps' between pos- 
itive and also negative values that correspond to open 
hyperbolic orbits. 

So an observer who assumes a pure Keplerian potential 
would deduce such an object moving on a hyperbolic orbit 
around the Sun if he lives in a MONDian universe. Whereas 
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Figure 10. Closeup of the semi-major-axis (SMA) between two 
peaks (from Fig. 9). One can see clearly that a takes positive and 
negative values, corresponding to closed and apparently hyper- 
bolic orbits. 
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Figure 12. The behaviour of the instantaneous best-fitting (pas- 
sive) eccentricity with time under the influence of a heliocentric 
DM-minihalo. For the initial conditions a Newtonian circular or- 
bit with 50000 AU radius was used. 
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Figure 11. The behaviour of the instantaneous best-fitting (pas- 
sive) semi-major-axis (SMA) and the aphelion distance with time 
under the influence of A_. For the initial conditions a Newtonian 
circular orbit with 50000 AU radius was used. 




1e+09 



2e+09 3e+09 
Time [Days] 



4e+09 



Figure 13. The behaviour of the instantaneous best-fitting (pas- 
sive) eccentricity with time under the influence of MOND. For 
the initial conditions a Newtonian circular orbit with 50000 AU 
radius was used. 



Dark-Matter and A_ would still produce closed orbits. 

Here one can see one of the most important differences 
in this investigation between MOND on the one hand and 
A_ and DM on the other hand. 

Fig. 12-14 show the passive eccentricities for a circular 
orbit in pure Keplerian case under the impact of the three 
modifications. 

The eccentricities in the DM and A_ cases take values 
betwen and 1 which corresponds to closed orbits. Not so in 
the MONDian case. Here the passive eccentricity has values 
in a range between and 5.5 corresponding to the results 
for the semi-major-distance. 



4-1.2 Active orbital parameters 

The calculation of the active orbital parameters was done 
for every orbit with initial distances between 5000 and 
50000 AU. It should be mentioned here that the effective 
aphelion and perihelion of the turn of a rosette-orbit do 
not define the axis of symmetry of the turn of the orbit. 
Therefore we used a reference-ellipse approximating one 
turn of the orbit. So the following orbital parameters are 



with respect to this referenc-ellipse. 

In order to calculate the active aphelion-migration of 
the generated rosette we took the coordinates of the aphe- 
lion of one turn and determined the angle of the aphelion 
with respect to the foregoing turn. The results for all three 
modifications are shown in Fig. 15 where we have fitted a 
curve to the data values which turned out to be highly de- 
pendent on the initial form of the orbit (e.g circular orbit, 
elliptical orbit, high or low eccentricities...) 

Very obvious is the coincidence of the curves of DM 
and A_. MOND behaves completely differently. 

A similar behaviour can be observed for the active ec- 
centrities of one turn that are shown in Fig. 16: DM and A_ 
seem to correspond very well. 



5 CONCLUSIONS 

We investigated the influence of three modifications of 
gravity. 

The first one was a Heliocentric Dark-Matter halo with 
a special density-profile, causing a distance-independent 
additional acceleration. 
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Figure 14. The behaviour of the instantaneous best-fitting (pas- 
sive) eccentricity with time under the influence of A_ . For the 
initial conditions a Newtonian circular orbit with 50000 AU ra- 
dius was used. 
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Figure 15. Comparison between the aphelion-migrations per day 
for all three modifications. 



The second modification was in Newton's second 
law, F = ma, called MOND, whereas the third and last 
modification was a negative time-dependent cosmological 
"constant", A_, that causes an additional attraction 
towards the centre of the considered system. 

The used Dark-Matter-halo had a special density pro- 
file ~ - that generates a constant acceleration towards the 
centre, which was indicated by the Pioneer-acceleration. 
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All the modifications caused strong deviations from 
Keplerian circular orbits. Circular orbits are transformed 
into elliptical ones experiencing in addition an aphelion- 
migration (rosette-orbit). 

The analysis of the passive and active orbital parame- 
ters has shown strong differences between MOND on the 
one side and A_ and Dark-Matter on the other one. The 
MONDian effect seems to dominate even for small distances 
of about 5000 AU whereas the other two modifications 
vanish at such distances. An explanation for this difference 
can be found in the non-linearity of MONDj and the fact 
that the contribution to the pure Newtonian acceleration is 
fairly larger than for A_ or DM. 

Although the Dark-Matter density-profile has been 
chosen such that it causes a constant acceleration like the 
Pioneer-acceleration neither this modification nor the two 
others investigated here can explain the Pioneer-anomaly. 
The distances where these modifications would have a 
measurable effect are far beyond Pioneer-like distances. 

However, every considered modification of gravity 
would cause effects that should be observable. In the future 
it would appear useful to do further investigations. 

Every modification of gravity has been considered at 
large distance of about 5000 — 50000 AU and seems to be 
negligible on scales of the inner Solar (« 50AU) System. 
A direct observation at distances of the order of the Oort 
cloud would be very difficult since every observation is 
limited by the possible resolution and sensitivity of the 
used telescope and furthermore by the observation time 
that would be necessary to observe deviations from the 
pure central mass Newtonian case at such large distances. 
However, if future-telescopes will be able to resolve e.g. Oort 
Cloud objects directly a comparison between the deduced 
orbital parameters with the passive (and active) parameters 
presented in this paper, would confirm or exclude the 
validity of the considered modifications of gravity. 

Another possibility besides the direct observation of 
long periodic Oort cloud objects could be an experiment on 
an Earth-bound satellite placed at one of the Lagrangean 
points. These points are quasi-stable points within the 
framework of the Three-Body-Problem where the net 
gravitational force of the bodies is almost zero. It is for 
this reason why these points would be appropriate to test 
the modifications due to A_ and a Dark-Matter halo. To 
test MOND these points would not be appropriate because 
the external Galactic field which has been neglected in this 
paper is much larger than an over the whole Solar System 
and so would not lead to any MONDian effects over the 
Solar System. Investigating the additional forces acting on 
an object at the Lagrangean points and comparing them 
with the forces due to A_ or DM would be a possible test 
of modified gravity. In addition also other modifications 
could be tested with such an experiment. Furthermore it 
could be useful to apply A_ on larger scales like binaries, 
disk-galaxies or dwarf-spheroidal galaxies. Also an investi- 
gation of the contribution of A_ to the perihelion-turn of 



Figure 16. Comparison between the active eccentricities for all 
three modifications. 



MOND deals with the absolute values of the acceleration 
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Mercury seems to be a worthwhile application. 

Whatever the investigation and development of modified 
gravitational theories will create, it seems to be obvious 
that the pure Newtonian Gravity due to central masses 
is not the whole story if the Pioneer-anomaly cannot be 
solved otherwise. 

One ansatz co uld be the one postulated by 

iBrownstein fc Moffatl (|2006T l. They postulated that the 
Pioneer-anomaly is caused by a distance dependence of 
Newtons gravitational constant G(r) on the basis of the 
metric-skew-tensor gravity (MSTG) and the scalar-tensor- 
vector gravity (STVG). It seems that this ansatz could be 
succesfull: "We have demonstrated that the STVG theory 
can explain the Pioneer anomalous acceleration data and 
still be co nsistent with the accur ate equivalence principle 
[...]" (JBrownstein fc Moffatl (|2006h . p. 3435). 
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